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\ ' , Abstract 
> 

, Applying a variational method to a Gaussian wave ansatz, we have derived a set of semi-classical 

' evolution equations for SU(2) lattice gauge fields, which take the classical form in the limit of a 

' vanishing width of the Gaussian wave packet. These equations are used to study the quantum effects 
on the classical evolutions of the lattice gauge fields. 
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1 Introduction 



Chaotic behavior in non-abehan gauge fields is first shown for the spatially homogenous potential config- 
urations jl] and later also for the spherically sysmmetric fields . Recently a real time evolution method 
is introduced to study general non-abelian gauge fields on a 3 dimensional lattice, where chaoticity is 
shown by obtaining a positive maximal Lyapunov exponent The maximal Lyapunov exponent is 
found to be related to the damping rate of gluons at rest, an infrared property of the quantized system. 
It can be argued that in the high temperature limit the properties associated with the infrared degrees 
of freedom in the gauge system are classical at the leading order. But it is still reasonable to ask how 
reliable these classical results are and what kind of effects the quantum corrections might bring about. 

A classical dynamical system can be studied in two ways. One is by tracing single particle tra- 
jectories using Newton's equation. The other is by studying the phase space density evolution based 
on Liouville's equation. The links between classical and quantum physics have been studied in both 
formalisms for a long time. One of these links is the Wigner function which is the counterpart of the 
phase space distribution in classical dynamics |^ . The driving force of its evolution can be separated into 
classical and quantum parts ^ . Unfortunately the variables in the lattice gauge theory are defined on a 
compact group manifold with a non-trivial geometry. The Wigner function on a compact position space 
is defined on a discrete momentum space rendering the classical limit rather subtle. Therefore we 
want to explore a more direct analogy to a classical trajectory in quantum physics, namely the evolution 
of a narrow wave packet. In order to resemble a classical trajectory as close as possible, the width of the 
wave packet in coordinate space and that in momentum space shall satisfy the uncertainty relation most 
efficiently, which suggests that we shall choose a Gaussian. 

In this paper, as our first attack on the problem of quantum effects on the classical evolution of 
a gauge system, we apply the variational method to Gaussian wave packets ^. The basic idea is to 
parametrize a wave packet, usually a Gaussian one, by a small number of parameters. Then the time 
dependent variational principle is used to derive a set of evolution equations for these parameters, which 
have a symplectic structure like Hamiltonian equations. If the parameters are chosen in such a way that 
some of them resemble the variables in the corresponding classical system, then the coupling between the 
semi-classical time evolution of these parameters and that of all others can be interpreted as quantum 
corrections to the classical evolution. This method has been used to treat the scattering problem of 
a particles [|9|, and recently it also has been applied to the many body problem of interacting fermion 
systems jiol, In the following we will first introduce the method. Then it will be applied to the pure 
SU(2) lattice gauge system. The semi-classical Hamiltonian and the equations of motion will be derived 
and their solutions will be studied numerically. 



2 Method 

We begin by considering a quantum system described by a Hamiltonian H . The dynamics of the system 
is contained in the Schrodinger equation 

indAm) = H\m)- (1) 

We solve this equation by the time-dependent variational method [|j as follows. A normalized trial 
wavefunction |$(q)), where q denotes a set of variational parameters, must satisfy the following variational 
principle 

5 dt{<^{q{t))\{ihdt~H)\<^{q{t)))^0. (2) 
Jtl 

Among these parameters some have definite meanings in the classical limit, such as the classical momenta 
and positions. Provided that enough parameters are chosen to span the whole Hilbert space, (2) can solve 
(1) exactly. Another case, when this is possible, is if we choose the parameters so clever that the exact 
time evolution of the initial state admixes only states from that part of the Hilbert space which is spanned 
by these trial states. But generally (2) is only an approximation to (1). Performing the variation in the 
usual way we can derive an effective Lagrangian 

L = qMq)-'H{q) (3) 
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^i(9) = ^(<^ll? (4) 



with 

2V dq.. 
and 

niq) =< > . (5) 

In the equations of motion, 

Aij{Q)Qj = 1^, (6) 
the antisymmetric coefficient, 

(7) 



q'=q 



occurs, which satisfies the Jacobi identity as well. 



OA- 

-j-^ + cyclic permutations = 0. (8) 
dQk 

Since Aij is antisymmetric, the expectation value of the energy, H, is conserved, 

— = qiAijqj = 0. (9) 

Other conservation laws related to symmetries of the quantum system can be derived accordingly. 

Suppose that G is a generator of an infinitesimal symmetry transformation related to a conserved 
charge, i.e. it commutes with the Hamilton operator 

[H,G]=0. (10) 

Denote the corresponding 'classical operator' by Go, whose action on the respective classical parameter 
xo is identical to that of the 'quantum operator' G on the variable x 

G:x^ fix) Go:xo^ f{xo). (11) 

We call our trial wave function covariant under the transformation G if it is invariant under the joint 
operation of G and Go, 

(G + Go)|$>=0. (12) 
Such wave functions lead to a conserved expectation value 

with 

e;=<$|G|$>. (14) 

To prove this, let us consider an infinitesimal symmetry transformation generated by Gq. From 
the variational principle (2), we have 

j dt (< 5<^\{ihdt -H)\<^> + < <^\{ihdt - H)\S^ >) = 0, (15) 

with 

>= iSaGol^ > 

< 6<^\ = -iSa < *|Go, (16) 
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where 5a is a small real parameter. 

The above integrand can be casted into the form 

-5a< ^IGol^* > -5a < $|Go|$ > -iSa < $|[Go,i/]|$ > . (17) 

Now using Go|<& >= — G|<i> > the commutator term vanishes, because G commutes with H. The rest 
can be splitted into a total time derivative, which we ignore in the variational integral, and into terms 
proportional with Sa. Finally we get 

'dt(5a (< $|Go|$ > + < $|Go|$ >) = 0. (18) 

Using again Go|$ >= — G|$ > and the fact that G = 0, we arrive at 

- J dt5a (< $|G|$ > + < $|G|$ > + < $|G|$ >) = 0, (19) 

which, for 5a being arbitrary, is equivalent to 

^<$|G|$>=0. (20) 

We intend to use this property of the time-dependent variational principle in the following in order to 
ensure that Gauss's Law is satisfied on the average during the time evolution of Yang-Mills systems. 

3 SU(2) Lattice Gauge Theory 
3.1 Variational Wavefunctions 

Now we apply this method to lattice gauge theory. Here we restrict ourselves to SU(2) theory for 
simplicity. The lattice gauge Hamiltonian is 

H = -{y: ^E^E^ + > -(1 - ^trUp) I , (21) 



(^^iw + ^ReE(l-itr[/,)^ 



where the first term corresponds to the electric and the second to the magnetic energy. In the classical 
limit the pure gauge theory is scale invariant. The coupling constant and the lattice spacing a can 
be scaled out entirely and we are left with no free parameter except the total energy (or temperature). 
In the semiclassical approach we want to consider the quantum corrections to classical trajectories so 
Planck's constant h comes in naturally. After making the following scaling transformation, 

h' = ag^H, 
t = t/a^ 

E = g^E, 

h = g^h, (22) 

the variables on the left become dimensionless numbers. In the following all quantities are understood 
to be the scaled ones unless otherwise stated explicitly. The equations all stay the same except that we 
shall put a = 1 and g = I- The scaled parameter h = 47ras, where as = g^Ti/An is the dimensionless 
coupling constant in the gauge theory, enters our formalism in two ways. First it appears in the variational 
principle, 

5 / {<i>\{indt - H)|$) = 0. (23) 
Jti 

Secondly it also enters the commutation relation, 

[Er,Al] = -zh5''%m. (24) 
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The expUcit appearance of couphng constant in our calculation fixes the scale of quantum effects relative 
to classical effects. 

To proceed we need to choose a trial wavefunction. In the strong coupling and zero temperature 
limit, the ground state wave function is an expansion of plaquette variables trUp It is separable for 
different plaquettes, but this choices is not very practical in our context. First we want the trial function 
to have the feature that the averages can be performed in closed analytical expressions. A plaquette 
separable wavefunction is more difficult in this aspect than a link separable one (see below). The second 
reason is that we are mainly interested in the high temperature limit where the system can be treated 
classically. In this limit the former is not a priori better justified than the latter. Moreover, since in the 
classical limit we deal with link variables rather than plaquette variables, a link separable wavefunction 
allows to make a more direct connection to classical dynamics. We choose the following Gaussian ansatz. 



(25) 

I rv^vi vz / 

The Uio and Eio are parameters of the wavefunction, which correspond to classical link variables and 
left electric fields, respectively. Physically they denote the center of the coordinates and momenta of the 
Gaussian wave function. The Ui's and J7;o's are group elements of SU(2), whereas the Eiq are elements 
of the su(2) algebra. The complex parameter b = v + iw controls the width, or more explicitely, is 
the inverse of the width. Ni is a normalization constant to be specified later. 

The form (25) is chosen so that it is covariant under the gauge transformations, i.e. invariant under 
the simultaneous gauge transformation on both the quantum variables Ui and the classical parameters 
Uio and Eiq, 

U'i=gLUig^^ U'lQ ^ gLUmg^^ and E'l^ = giEmgl'^ , (26) 

where gL{R) is a gauge transformation at the site left(right) to the link I. This leads to the desirable 
result that Gauss's law is conserved on average. 

We note here that color confinement requires that all physical states shall be color singlet, while 
the trial wavefunction we use here is not gauge invariant. Is this a contradiction? The answer is yes in the 
sense that our choice does not correspond directly to a physical state. But here we have two arguments to 
support our choice. First we note that a covariant discription is as close as we can get in the classical limit 
where we are always using color octet variables Ui and Ei to define a gauge configuration. Had we started 
from a color singlet trial wavefunction and wanted to derive classical equations of motion for averages of 
these variables, we would always get vanishing expectation values because of the Wigner-Eckart theorem. 
The second argument is that given a gauge covariant wavefunction we can easily project out its singlet 
part. So the knowledge of the former can help to understand the properties of a gauge invariant state. 

In the continuum limit, this wavefunction formally coincides with a Gaussian in the vector potential 

Aix), 



$ [{A{x)] = exp 



(27) 

The wavefunction (25) is normalized on each link with 

N{v) ^2t:^Ii{2v)/v, (28) 

where /„ denotes the modified Bessel function. Wavefunctions with different parameters are not orthog- 
onal to each other. They form a set of overcomplete basis states of the Hilbert space for each link, with 
the completeness relation 

dUod^'EoHEo, C/o, bMEo, f/o, b) = 1, (29) 

for arbitrary b. This permits us to deal with the complex width parameter b in two ways. We could 
define a measure /i(5) which satisfies 

/ rfM(^) = 1, (30) 
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so that we formally have complete closure relation over all parameters 

dUod^EodfxibMEo, Uo, b)^{Eo, Uq, b) = 1. (31) 



The exact form of /i(5) is not known, but could be obtained with numerical techniques. On the other 
hand we can treat b as an additional variational parameter. We will adopt this point of view here. 
The relation between our classical parameters and the averaged quantities are obtained as 

(S-) = fiv)ES, 
(i^S) = f{v)E%o, 

{U) = f{v)Uo, (32) 
where -E^r" = Ut^U~^E^ are the so called right electric fields and the function, 
.M_'^ 9lnN{v) _ h{2v) / {v /2){l - /Q) v^l 

•^^ ^~2 dv ~ h{2v) \ l-3/(4u) u»l ^"^"^^ 

is plotted in Fig.l (solid lino). Due to the combined effect of the curvature of the SU(2) group manifold 
and the finite width of the wave packet, the average of E and U are not just their corresponding classical 
parameters as we would expect for a Euclidean parameter space. In the small width limit, i.e. v ^ oo, 
f{v) approaches unity, a manifestation that locally the group manifold has Euclidean metric. From 
this we already see that the non-trivial geometry plays an important role here. We may also study the 
quadratic fluctuations of U and E. We calculate 



{E^) = Y.iE'^E^ = h'^^iv' + W-) El 



2 (34) 

-1 



We shall later see that the ensemble average of v is propotional to and ensemble average of w is zero 
for V ^ oo. So in the classical limit, i.e. h ^ 0, {E"^) goes to the classical result Eq. The variance of 
electric fields is 



which goes to zero like + ^ in the classical limit. The fluctuation of the vector potential A is 

{AAf ^ 8(1 - ^tr([/C/-i)) = 8(1 - /). (36) 

It vanishes in the classical limit like h. We can check the uncertainty relation AEAA > 3{h/2), which 
becomes an equality in the vacuum, i.e. T = or E^o = 0. As we already said. Gauss's law is satisfied on 
average, but it is not conserved at the operator level. The magnitude of fluctuating color charge at each 
lattice site is measured by 

{AGf = ((G^G-^) - (G")2) , (37) 

a 

where the second term vanishes. Inserting the expression of the generator of time independent gauge 
transformation, 

3 

G'' = ^(£;"(n)-ES(n-e,)), (38) 

2 = 1 

where the summation runs over the three positive directions at one site, and n denotes the position of 
the site and ej is the unit vector in ith direction, we obtain 

{AGr=j2i^E)h (39) 

i=l 

where the summation now runs over all the links joining at one site. Each link contributes its electric 

field fluctuation to the local charge fluctuation. The correlations between different links do not appear 
because we have chosen a link separable wavefunction. 
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3.2 Effective Hamiltonian 



The complete expression for the average total energy is 



Hsc = He + Ho + Hm = {H) 

^1 



E 

I 



2vi 
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E 



-JptvUpo 



3.fi 
16vi 



ivf+wf) 



+ 



(40) 



where fp is the product of fi = of the four links that circumvent the plaquette p. This semi-classical 
Hamiltonian consists of three parts. First we have the modified electric energy {He) and magnetic energy 
{Hjn), which result in a coupHng between the Ui, Ei and the width parameters vi and Wi, and thus 
modify the classical evolution of Ui and Ei. These modifications contain both genuine quantum effects 
and lattice artifacts, which are difficult to separate. Besides the modified classical electric and magnetic 
energy there is a zero point energy (Hq) which is a consequence of the uncertainty relation. It increases 
linearly with v in the limit of large v. 

We can study thermodynamical properties of the system described by Hsc- Consider the system 
at equilibrium with temperature T. The integration measures of E and U are the same as in the classical 
case. The integration over E is plagued by Gauss's law which eliminates one third of the degrees of 
freedom. The average magnetic energy is complicated by the interactions between different links. They 
make it very difficult to obtain the exact solution analytically. In the rest of this subsection, let us 
simplify ourselves to consider the Hamiltonian of a single link and a single plaquette. In this free link 
and plaquette limit the partition function can be factorized into three parts, 



Zf,Zm.Zl 



with 



0) 



T 



(41) 



3/2 



Zo 



7V(2/p/T)exp(-4/T), 
eM-Ho/T), 



where the subscripts e, m, denote contributions from electric, magnetic and zero-point energy, respec- 
tively. N{x) is the normalization function defined in (28). From the partition function we can easily 
obtain the ensemble averaged energy 



{{He)) = 3T/2, 
{{Hm)) = 4(1 - /p/(2/p/T)), 
{{Ho)) = Ho, 



(42) 



where ((• • •)) denotes ensemble average. We see that the electric energy is the same as that in the classical 
limit. In the classical Hmit fp 1 so {{Hm)) = 4(1 - f{2/T)). When T < 2, {{Hm)) takes the correct 
classical limit 3T/2 as demanded by the equipartition theorem. At larger T the energy temperature 
relation is distorted by the compactness of the gauge group. The free energy F = —TlnZ is obtained as 



XT 

F{v) = —\n 



^ v'T 



+ 



16 



fv + A-TlnN ' ^-^^ 



(43) 



To see how the free energy looks like, we plot this function at T = 0.5 and at T = in Fig2. The two 

plots are very similar in shape, and both have a first order phase transition. At small h the free energy 
has two minima, one of which is at = 0. The second minimum can be obtained in the limit of large v 
by minimizing the free energy, and the result is 



V = 2(16 - ITfl'^l'h. 



(44) 
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When h increases, this second minimum rises until it fails to be a minimum at a critical value of h, where 
a first order transition occurs. The critical value increases with temperature. From the plots we see for 
T = it occurs at h ^ 1.73 (the short dashed line in Fig. 2b) and h ^ 3 for T = 0.5 (the long dashed line 
in Fig. 2a). This indicates the underlying potential has a non-trival structure. Vacuum energy is obtained 
from the free energy by setting T = 0, 

Fo = —fv + i{l~n. (45) 

At large h, we only have one minimum v = 0. The wavefunction is so broad that it covers the whole 
group manifold and we are in quantum physics region. The corresponding energy is 4 (or i/g^a in original 
variables). This is coincident with the leading order result of a perturbative calculation in the strong 
coupling limit [p^ . In the other limit, ?i ^ 0, the wave packet can be arbitrarily narrow, i.e. v —> oo, so 
it can probe the lowest point of the potential well, and we are in the classical regime. The energy density 
in this limit vanishes. 



3.3 Equations of Motion 

To obtain the semi-classical equations of motion we calculate the expectation value of the time derivative: 

(ihdt) = J2^~nfwi + iftriEioUioUi-Q^)). (46) 
I 

Inserting this expression as well as the one for (H) in (2) and performing the variation, we obtain the 
equations of motion for each link (the subscript is suppressed and fi stands for f{vi)): 



Vl 



Wl 



^(^--t) , (47) 



where /; is the derivative of /; with respect to w/, the sum over p{l) runs over all four plaquettes sharing 
the link I and Up here is the ordered plaquette variable starting with link [//. We note that in the limit 
of small width or large v the equations assume the correct classical form, 

E- = z^tr(T'^C/p), 

Ui = lEiUi. (48) 

Hence the equations (47) can be used to studied the quantum corrections to the classical evolution (48). 
These equations can also be studied in the limit where v is fixed and w = 0, where the equations take 
the form 

p(i) 

Ui = 7^{l-^)EiUi. (49) 
2/; 2u; 

The fixed width parameter v is to be treated as a variational parameter. 
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3.4 Classical Limit 



In order to bo ablo to define a trajectory in phase space the wave packet must be narrow compared to the 
total available phase space volume. The restriction comes both from the link variables U (coordinates) 
and electric field variables E (momenta). To have a quantitative measure of how localized the wave 
packet is, we define a quantity a{v) which is the ratio of the volume of the wave packet and the total 
group volume, 2tt^ for SU(2), i.e.: 

a{v) = {l-ff/'^\ ^3,3/2 """^ (50) 



(to) 



oo. 



The function a(v) is shown in Fig.l (dotted line). If we request a{v) < 0.1, we need v > 5 and h < 1 
for T « 1. The largest distance in momentum space is controlled by excitation energy as d^; ~ ^/E or 
ds ^ T. We require AE <C cZb- In the classical limit {v ^ oo) this requires h <^T, which in the unsealed 
variables is 

A = ^ « 1. (51) 

where kmax is the largest wave number that is possible to appear on the lattice. Physially this means 
that we shall only keep the long wavelength modes on the lattice which is consistent to our intuition. 



4 Results and Discussions 

With the above semi-classical equations of motion we can follow a trajectory in the space of gauge fields 

for every given initial condition. To make the results comparable to those of our calculations for the 
classical gauge theory, we specify the initial gauge configuration in the same way as in ref. [3] . We initially 
choose all electric field variables to be zero to satisfy Gauss's law. The magnetic sector is initialized 
by chosing the link variables Ui = cos(p/2) — m • Tsin(p/2) randomly and the total energy is varied by 
selecting 6 which limits the range of the parameter p to (0,27r5). In addtion here we specify h and the 
width parameters. Initially we choose identically on each link v = 20 and w = 0. The results are given 
in Table 1. 

To present our result we need to specify our controllable parameters. One of these is h = 47ras, 
which is a manifestation of quantum effects. Another should be chosen to give a measure of the excitation 
energy of the system. A direct comparison between classical and quantal systems with the same energy 
is unphysical because of the appearance of a zero-point energy in the latter. A convenient comparable 
parameter is the physical temperature of an excited system. The problem now is how to measure the 
temperature on the lattice. For a weakly interacting system, we can extract the temperature as follows. 
The energy distribution for a small subsystem is measured, which factorizes into a phase space prefactor 
and an exponential ex.p{—E/T). Dividing out the phase space factor we can simply identify the temper- 
ature as the rate of exponential slope of the energy distribution. But here this method cannot be used 
directly because the lattice gauge system is a constrained and strongly coupled system. This problem 
will be discussed in detail elsewhere. Here we make the observation that the relative difference between 
the electric part of the semi-classical Hamiltonian and that of the classical one is l/(2t;) which is very 
small in the limit of large v. As a consequence, the value of {{He}), the average electric energy per link, 
at a certain temperature T of the classical system and that of the semi-classical system in the limit of 
large v will be similar. Neglecting the small difference, we can use it as a measure of the excitation level 
of the system, and define parameter Tg = {{He)) as the "electric temperature" of the gauge field. 

We have measured width parameters on the lattice. The results are shown in Fig. 3 for h = 0.3 
and T(. = 0.9, and in Table. 1. In Fig. 3a and 3b, the time evolution of the average values of v and w are 
shown. The average is taken over the lattice, which is expected to coincide with the ensemble average 
if the lattice is siifficiently large. We find that these values are quickly damped into their equilibrium 
positions, which are independent of their initial value. The final distributions of v and w over the lattice 
are shown in Fig. 3c and 3d respectively. The imaginary part of the width parameter w obviously has 
zero mean value. The nonzero expectation value of v depends mainly on h as we can see from column 4 
of Table 1 . We have obtained the least error fit 

{{v)) « 7/n. (52) 
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The T,. dependence is difficult to see because T,. does not change much between runs. The fi dependence 
roughly agrees with (44), which is found for the system without constraints and couplings by the standard 
thermodynamical method. 

We have also measured the electric energy distribution which is shown in Fig. 4 for the same initial 
condition as in Fig. 3. The histogram is from the simulation. The solid curve is a fit, 



apart from a normalization constant. It has correct mean value {{He)) = 0.9 and looks very much like a 
Boltzmann distribution. We note that for a constrainted system it is possible for the exponential in the 
energy distribution to have the form exp{—(3E/T) where /3 can be different from 1. So here we do not 
conclude the system temperature is 0.6. 

Now we come to the question which originally motivated this study, namely, whether the quantum 
effects suppress or enhance the chaotic behavior of a pure gauge system. For this purpose we shall 
measure the maximal Lyapunov exponents A in our semi-classical formalism and compare them with the 
ones we found for the classical lattice gauge theory at the same Tg. The logrithmic scale of the divegence 
of trajectories is shown by the solid line in Fig. 5, where h = 0.3 and = 0.9. For comparison, we also 
show in the same plot the divegence of two corresponding classical trajectories, with same Te, by the 
dotted hue. The distance measure used here is the same as that in ref.[3], i.e., the integrated absolute 
value of the local difference in magnetic energy. At early times, while the classical evolution smoothly 
increases the distance between the two trajectories, there is a sudden jump of the distance to a value 
independent of Dq in the semi-classical evolution, which originates from the couplings of Ui and Ei to 
the width parameters vi and wi . The latter oscillate and are damped strongly in the beginning as we see 
in Fig. 3. Also at early times the fluctuations are large on the logarithmical scale compared to the purely 
classical evolution. These fluctuations are quantum effects. At larger times the exponential increase of 
the distance emerges from the fluctuating background. We see the rate of exponential growth is larger 
than that of the classical one. The distance finally reaches a maximum value. The maximal Lyapunov 
exponent is extracted from the linear region in the plot. Due to the large fluctuations this method is 
less accurate here than it is in the classical case. The numbers obtained for A have errors of the order of 
±0.02. We performed three nms for the system with the same h and Te and find that the statistical error 
is dominated by our measuring error. The results for different initial conditions arc shown in column 3 
of Table. 1. At ?i = 0.1 the result deviates very little from the classical result Ac = Te/3. The correction 
to A is represented by d\ (column 4) which is defined as X{Te, Ti) = Ac(Te)(l + SX). SX depends both on h 
and Tg. The relation is not easily extracted from the data due to the lack of accuracy, but we note that 
the sign of 6 is always positive. This suggests that quantum corrections enhance the chaoticity of the 
classical gauge system. The physical implication is that including the quantum corrections, the initial 
state of the gauge field will tliermalize more rapidly. 

We have also measured the time dependence of the distance between two trajectories originating 
from the same point but one evolved with the classical equations and the other with the variational 
quantum equations. This distance also diverges exponentially, because as soon as the quantum corrections 
generate deviation between the trajectories the difference will diverge with time because of the chaoticity 
of our system. This leads to the question whether the trajectory of the centroid of the wave packet has 
any physical meaning. If we want to describe the quantum corrections to a particular classical trajectory 
exactly, then the method of Gaussian variational states obviously fails. However, if we are only interested 
in the ensemble averaged properties of the system, e.g. the rate of the increase of the coarse grained 
phase space volume, the approach remains useful. This is not too serious a drawback because even in 
the classical limit we cannot numerically follow the exact evolution of a particular trajectory for a long 
time, and we have to limit ourselves to the study of average properties of the system. 

In conclusion, starting from a Gaussian ansatz we have derived a set of semi-classical equations of 
motion for the lattice SU(2) gauge field via a variational method. These equations were used to study 
the quantum corrections to the classical evolution. We find the quantum effects enhance the chaoticity 
of the gauge field. 
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Table Captions 

Table. 1 Results of numerical simulations. 

Figure Captions 

Fig.l The solid line shows function f{v) and the dotted line shows wave packet volume ratio C{v) as a 
function of v . 

Fig. 2 Dependence of the free energy on v for different values of ?i at T = 0.5 (a) and T = (b). The sloid 
line is for h = 0.1, the dotted line is for 1, the short dashed line is for 1.732 and the long dashed 
line is for 7i = 3. v ^ oo corresponds to the classical limit of an infinitely narrow wavepacket. 

Fig. 3 Width properties for Tg = 0.9 and H = 0.3: The equihbrium distribution of v (a) and w (b) and 
time evolution of the mean values of v (c) and w (d). 

Fig. 4 Final electric energy distribution for Tg = 0.9 and ?i = 0.3. 

Fig. 5 The solid line shows the divergence of the centroids for the two neighbouring Gaussian field con- 
figurations for Te = 0.9 and h = 0.3. The corresponding divergence for the two classical field 
configurations with the same Tg is shown by the dotted line. 
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Table. 1 
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0.24 


0.12 


23 


2.48 


0.3 


0.73 


0.3 


0.22 


22 


2.92 


0.3 


0.85 


0.34 
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21 
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0.85 


0.36 


0.27 
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3.25 
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0.41 


0.37 
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3.74 
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0.19 


0.16 
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0.41 


0.54 
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4.40 
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